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We consider the possible existence of self-similar solutions containing black holes in a Friedmann 
background with a stiff fluid or a scalar fleld. We carefully study the relationship between the 
self-similar equations in these two cases and emphasize the crucial role of the similarity horizon. 
We show that there is no self-similar black hole solution surrounded by an exact or asymptotically 
flat Friedmann background containing a massless scalar fleld. This result also applies for a scalar 
fleld with a potential, providing the universe is decelerating. However, if there is a potential and 
the universe is accelerating (as in the quintessence scenario), the result only applies for an exact 
Friedmann background. This extends the result previously found in the stiff fluid case and strongly 
suggests that accretion onto primordial black holes is ineffective even during scalar field domination. 
It also contradicts recent claims that such black holes can grow appreciably by accreting quintessence. 
Appreciable growth might be possible with very special matter fields but this requires ad hoc and 
probably unphysical conditions. 
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I. INTRODUCTION 

Scalar fields are one of the key ingredients in modern 
cosmology. They play an important role during inflation 
and certain phase transitions; they feature in the preheat- 
ing and quintessence scenarios; and they are pervasive in 
all sorts of alternative theories (eg. string theory and 
scalar-tensor theory) which are likely to be relevant in 
the high curvature phase of the early Universe. In such 
scenarios, it is often assumed that scalar fields dominated 
the energy of the Universe at some stage of its evolution. 
This could have profound implications for the formation 
and evolution of any primordial black holes (PBHs) that 
may have formed in the early Universe [J . It would also 
modify the constraints on the number of PBHs provided 
by various observations 0, Q . 

A particularly interesting example of this arises if 
the dark energy which generates the cosmic accelera- 
tion observed at the present epoch is associated with a 
quintessence field. In this context, it has recently been 
claimed jj, yf that PBHs could grow enough through 
accreting quintessence to provide seeds for the sort of 
supermassive black holes thought to reside in galactic 
nuclei 6]. This claim is based on a simple Newtonian 
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analysis, in which a PBH with size comparable to the 
cosmological particle horizon appears to grow in a self- 
similar manner [3, i.e. its size is always the same fraction 
of the size of the particle horizon. This suggests that 
PBHs may accrete surrounding mass very effectively if 
their initial size is fine-tuned. However, this analysis ne- 
glects the cosmological expansion i3,ip|. More recently, a 
general relativistic analysis of this problem has been ap- 
plied for a black hole accreting dark energy and phantom 
energy [iflll3)ll3> ^^ well as a ghost condensate [l3l[l4l| . 
However, the cosmological expansion is still neglected in 
these analyses. 

The possibility of spherically symmetric self-similar 
PBH solutions was first studied more than 30 years ago 
for a fluid with a radiation equation of state {p = e/3) Q]- 
This study was also prompted by the Newtonian analy- 
sis and demonstrated analytically that there is no self- 
similar solution which contains a black hole attached to 
an exact Friedmann background via a sonic point (i.e. in 
which the black hole forms by purely causal processes). 
The Newtonian analysis is therefore definitely mislead- 
ing in this case. There are self-similar solutions which 
are asymptotically Friedmann at large distances from the 
black hole. However, these correspond to special initial 
conditions, in which matter is effectively thrown into the 
black hole at every distance; they do not contain a sonic 
point because they are supersonic everywhere. 

Similar results were subsequently proved for all perfect 
fluid systems with equation of state of the form p = ke 
with < fc < 1 |l^ 01 . Indeed it was shown that 



the only physical self-similar solution which can be at- 
tached to an external Friedmann solution via a sonic 
point is Friedmann itself; the other solutions either en- 
ter a negative mass regime or encounter another sonic 
point where the pressure gradient diverges 16]. Later 
numerical relativistic calculations have confirmed that 
PBHs cannot accrete much for these perfect fluid sys- 
tems ,18;, 19, 20, 21, 22, 23, 24]. It has also been shown 
that the asymptotically Friedmann self-similar solutions 
with black holes are not strictly Friedmann at infinity, 
because they exhibit a solid angle deficit which could in 
principle show up in the angular diameter test [25j. 

The limiting values of k require special consideration. 
In the fc = case, there is no sonic point, so the unique- 
ness is trivial; again there are asymptotically Friedmann 
self-similar solutions with PBHs but only for special ini- 
tial conditions Q. The stiff fluid case (/c = 1) is more 
complicated and has an interesting history. It was orig- 
inally claimed by Lin et al. ^26] that there could exist a 
self-similar black hole solution in this case but subsequent 
calculations by Bicknell and Henriksen z7] showed that 
this solution does not contain a black hole after all be- 
cause the alleged event horizon is timelike. The only way 
to avoid this conclusion would be if the fluid turned into 
null dust at this timelike surface, but this seems physi- 
cally implausible. It therefore appears that there is no 
self-similar solution containing a black hole in an exact 
Friedmann background for any value of k in the range 
[0,1]. 

In view of the prevalence of scalar fields in the early 
Universe, it is clearly important to know whether the 
non-existence of a self-similar solution extends to the 
scalar field case. One might expect this from the stiff fluid 
analysis, since a scalar field is equivalent to a stiff fluid 
provided the gradient is everywhere timelike |29 ]. This 
is supported by recent numerical studies of the growth 
of PBHs in a scalar field system, which - for a variety 
of initial conditions - give no evidence for self-similar 
growth ,2Qi 30]. However, it is very difficult to deter- 
mine through numerical simulations whether self-similar 
growth is impossible for any specific initial data. One 
therefore needs an analytic proof that there is no self- 
similar black hole solution if the cosmological expansion 
is taken into account. This paper presents such a proof 
and also provides the basis for studying the growth of 
black holes in an expanding universe containing more 
general matter fields. 

The plan of this paper is as follows. In Section II, we 
review the stiff fluid case, with particular emphasis on the 
earlier work of Lin et al. |26| and Bicknell and Henrik- 
sen J23. In Section HI, we consider the massless scalar 
field case, using Brady's formulation of the self-similar 
equations as an autonomous system [31i] and making a 
detailed comparison with the stiff fluid case. We prove 
that there is no self-similar solution unless the behaviour 
of the matter is very contrived. In Section IV we extend 
the analysis to the case of a scalar field with an exponen- 
tial potential (as applies in the quintessence scenario) and 



show that the proof applies in this case also. We sum- 
marize our conclusions in Section V. Various technical 
issues are relegated to appendices, including a compar- 
ison with previous literature, a discussion of the global 
characteristics of these solutions, and an analysis of the 
behaviour of solutions at the similarity horizon. We use 
geometrised units with c = G = 1 throughout this paper. 



II. STIFF FLUID CASE 

A. Bicknell and Henriksen's formulation 

First, we show how a stiff fluid description of a scalar 
field can give insight into the fluid-field correspondence 
for self-similar solutions [221 ■ ^ ^^^^ fluid is a perfect 
fluid with the equation of state p = e and stress-energy 
tensor 



T 



ah 



e(2ii°?i'' + .9°") 



(2.1) 



If M° is vorticity-free, one can show that this matter fleld 
is equivalent to a massless scalar field (f), for which the 
stress-energy tensor is 



Tn. 



■^9ab<f>,c(f 



(2.2) 



and the gradient 4'^a is timelike. The associated energy 
density and 4-velocity are 



= ±- 



V-^^ 



(2.3) 



where the sign in the second equation is chosen so that 
u"' is future-directed. The line element in a spherically 
symmetric spacetime can be written as 



-e2-(*.«)di2_^e2^(*'^)di?2-Hr2(i,i?)(d0' 



2 ' sin2 



ds' 

(2.4) 
where R and r are the comoving and area radial coordi- 
nates, respectively. 

Let us now assume self-similarity (i.e. the homothetic 
condition), which implies that we can write 



iy = u{X),X = X{X),S{X)^-, 



cj)^h{X) -Klnjtl, 



(2.5) 

where X = \n{R/\t\). The constant k is arbitrary but we 
will need to assume n — Xj \fV2/n below in order to allow 
the Friedmann solution. If we adopt comoving coordi- 
nates, which is always possible providing the 4-velocity 
is timelike, then T^^ is diagonal. In this case, hiX) must 
be constant and so^ depends only on t. 

Following Ref. [l£j, we now introduce the following 
variables: 



y2=e2A-2.+2X E = %^eR^, A^ 



Aner^ 



(2.6) 



where V is the speed of the fluid flow relative to the 
similarity surface X=const, m is the Misner-Sharp mass. 



£^ is a dimensionless measure of the energy density and 
A is one third of the ratio of the average density within 
r to the local density there. Energy conservation yields 
the following integrals: 
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axE-^S- 



(2.7) 



where a^ and a\ are constants of integration. Eqs. H2.3|l . 
(E!3) and (E21) give 



32"^ = 47r«;2s-ie2^, 



(2.8) 



so fly is identified as Attk^ and is invariant under rescaling 
of t and R. If we rescale R appropriately, we can also set 
ax to be ^ttk^. Eqs. (|TH|l and (|T7jl then imply [23| 



1/^ 



(2.9) 



so the Einstein field equations reduce to the following 
autonomous system j27.. .32. ] : 



{V')' = 2{l-A)V\ 
A' = \{\ + A){l-iA) 



-2AV 



E' = 2 



V' 



,A -Attk^ 
A ~ Attk'^ 



1 



E, 



(2.10) 

(2.11) 
(2.12) 



i-y2 

together with a constraint equation 

V^{l-Af~{l+Af+16TrK^V'^{E~^\V\-A) = 0. (2.13) 
Here a prime denotes a derivative with respect to X. 



B. V^ 



1 surface 



Once we have determined initial values satisfying 
Eq. f^TT^ . Eqs. (|TTU|) and (|CTT) const itute a two- 
dimensional autonomous system and Eq. (|2.12|l can be 
evolved freely. However, Eq. (|2.11(l shows that the condi- 
tion V"^ = 1 corresponds to a singular point in the above 
system of differential equations, so the usual uniqueness 
of solutions may not apply. This is because one may 
have a sonic point where V^ — 1 and, in the stiff fluid 
case, this coincides with a null surface; we term this a 
"similarity horizon" . 

In order to understand the significance of the similar- 
ity horizon, it is useful to draw the Penrose diagram for 
the self-similar black hole spacetime. This is shown in 
Fig. m which is also equivalent to Figure 5 of 33] . It can 
be obtained by combining the Penrose diagrams for a 
decelerating fiat Friedmann universe, which is discussed 
in Appendices IbI and O and that for a naked singular- 
ity |34| . The spacetime is asymptotically flat Friedmann 
with a spacelike big bang singularity outside the particle 
horizon. It also has a null naked singularity, a spacelike 
black hole singularity and a black hole event horizon. The 



black hole singularity 




.•0 



initial singularity 



FIG. 1: The causal structure of self-similar black holes in a 
decelerating universe. The trajectories of similarity surfaces 
are shown. There is a spacelike big bang singularity, a null 
naked singularity and a spacelike black hole singularity. The 
particle horizon (PH) and black hole event horizon (BEH) are 
both similarity horizons. 



similarity surfaces are spacelike outside the particle hori- 
zon and inside the black hole event horizon, timelike in 
between and null on the horizons themselves. This shows 
that the particle horizon and black hole event horizon 
both correspond to similarity horizons. 

We note that the induced metric on the X = const 
surface can be written as 



ds^ ^-e'^''{l-V^)dt^ + r^{d9^ + sin^ 



'■)■ (2.14) 



If this surface is null, it can be identified with the black 
hole event horizon or the cosnrological particle horizon, 
both corresponding to the condition V^ = 1 providing 
ly remains finite. If the gradient of the energy density is 
finite at the similarity horizon, Eq. (|2.12|l also requires 
A — Attk'^. This coincidence between singular points of 
the ordinary differential equations (ODEs), i.e. sonic 
points, and similarity horizons only applies for a stiff 
fluid. Sonic points are timelike surfaces for perfect fluids 
with the equation of state p — ke {Q < k < 1). 

If one wants a solution which represents a black hole 
surrounded by the exact Friedmann solution, we first 
note that the velocity gradient must be positive on both 
sides of the matching surface. (This is a consequence of 
the fact that the Friedmann sonic point is a node, as dis- 
cussed in Section lill El l Therefore V'^{X) must have the 
behaviour shown in Fig. [3 As X decreases from infin- 
ity, V'^ first decreases, goes below 1, reaches a minimum, 
then increases and reaches 1 again. The first point where 
V"^ crosses 1 is the particle horizon, X = Xph. One 
might expect the second point where V^ crosses 1 to be 
the black hole event horizon, X — Xbch, and indeed Lin 
et al. ^M\ make this interpretation. However, Bicknell 
and Henriksen ^27] showed that this is not the case: the 
similarity surface is timelike rather than null because v 
diverges there. More precisely, E goes to 0, which shows 
that the gradient of the scalar field is null, and the factor 
e2i'(l — V^) in Eq. (|2.14|) tends to a finite positive limit, 
which shows that the second crossing surface is timelike. 



This reveals the hmitation of the comoving-coordinate 
framework for scalar field systems. The Lin et al. and 
Bicknell and Henriksen analyses are discussed in more 
detail in Appendix IXI 




FIG. 2: The required behaviour of the velocity function 
V^{X) for a black hole solution which is attached to an exact 
Friedmann background. 

Because Bicknell and Henriksen could only find a so- 
lution in which V'^ has the required behaviour up to the 
point where V^ — 1, they inferred that no black hole 
solution is possible if the fluid has the same equation of 
state everywhere. However, they did find that a black 
hole solution is possible if the stiff fluid turns into null 
dust at the surface where V^ = 1. This situation is obvi- 
ously contrived, so in this paper we seek a more physical 
extension with a massless scalar field. 



III. MASSLESS SCALAR FIELD CASE 

A. Brady's formulation 

Let us now return to the original scalar field de- 
scription, in which the stress-energy tensor is given by 
Eq. (|2.2|) . If we adopt (non-comoving) Bondi coordinates 
(v,r), where v is advanced time, the metric can be writ- 
ten as 



ds^ = -ggdv"^ + 2gdvdr + r'^{d9'^ + sin^ 



'■), (3.1) 



where g and g are functions of v and r. This applies 
whether the gradient of the scalar field is timelike, null 
or spacelike. In this coordinate system, for a self-similar 
solution we can take the independent variable to be ^ = 
ln(r/|u|) and then define the following variables: 



y(0 



z{£,) = g ^x, 



g 2m 

g ^ 

cj,^m-n\n\vl 7(0 = MO, (3.2) 



where x = r/v = ±e^ has the same sign as v and a dot de- 
notes a derivative with respect to ^. Note that depends 
on ^ in this case, because we are not using comoving co- 
ordinates. The quantity z will play an analogous role to 
the quantity V in the stiff fiuid picture. The quantity 7 
is related to the gradient of the scalar field; the field itself 
never appears explicitly in the equations. 

With the self-similarity ansatz, the Einstein equations 
and equations of motion for the scalar field reduce to the 
following ODEs: 



(Ing)- = 



[ghy +gh 



^2 
47r/i , 

g - g, 

811(11 + K)[x(h + K) 

2x(h + h + 



gh], 



We then have an autonomous system 13 1|: 

z = z(2-y-i), 
y = 1 - (4^ + l)y, 
(l-2z)7 = 2Kz-7(y"^ - 2z), 



with the constraint equation 

y[{l + Attk"^) - 47r(7 + Kf{l - 2z) 



1. 



(3.3) 
(3.4) 

(3.5) 
(3.6) 



(3.7) 
(3.8) 
(3.9) 



(3.10) 



Note from Eq. (|3.2() that y < 1 if the Misner-Sharp mass 
is positive. This system was investigated in the context 
of cosmic censorship in l3,'l and, for k = 0, one has the 
exact Roberts solution i3S|. Eq. (|3.1U|) gives two values 
of 7 for each (y, z): 



-K± 




(3.11) 



with the positive (negative) branch corresponding to pos- 
itive (negative) values of 7-I-K. The structure of this sys- 
tem of equations is the same as in the stiff fluid descrip- 
tion. Since Eq. H3.11|l gives 7 in terms of y and z, once we 
have determined initial values satisfying the constraint 
(|TTn|) . Eqs. lEiZIl and ^^ together with Eq. (|XTT|) form 
a two-dimensional dynamical system and Eq. H3.9|l is 
evolved freely. 

Special physical significance is associated with the gra- 
dient of the scalar field being null. This is because the 
scalar field can no longer be identified with a stiff fiuid 
once the gradient has turned from timelike to spacelike. 
Since this gradient is 



^''''<i>,c = -7y[2z(7 + k) - 7]r 



(3.12) 



it can become null under three conditions: y = 0, 7 = 
and 22(7 + k) = 7. From Eq. (|3.1U|) . the condition 7 = 
gives 



1 



V 



SiTK-^y 



(3.13) 
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FIG. 3: The solutions in the yz plane for Attk^ — 1/3. The shaded regions are prohibited for a real scalar field. The expanding 
flat Friedmann solution is plotted as a thick solid line, the curve 7 = as a thin solid line, and the curve 22(7 + k) = 7 as 
a thin dashed line. The gradient of the scalar field can change from timelike to spacelike on the last two curves and also on 
y = 0. Note that the parts of the 22(7 + k) = 7 and 7 = curves in 2 > 1/2 correspond to the null conditions for different 
solution branches but are shown together for convenience. The flat Friedmann solution has a cosmological apparent horizon 
at {y = 0, z = 1/2), a particle horizon at {y = 3/4, 2 = 1/2) and ends a,t {y = l,z — 0). Two numerical solutions within the 
particle horizon are shown by the thick dotted lines. The upper one is a positive branch and goes directly to the negative mass 
region {y > 1). The lower one reaches y = 3/4, where it is converted from the positive branch to the negative branch (shown by 
the thick dashed-dotted line) before going to the negative mass region. The Bicknell-Henriksen null-dust black hole solution, 
which connects to the lower numerical solution, is indicated by a thin dashed-dotted line, but this is never realized by a scalar 
fleld. The black hole has an event horizon where the solution intersects 2 = 1/2 and an apparent horizon at {y — 0,z — 00). 



which always corresponds to the positive branch of 
Eq. H3.11|l . From Eq. H,'110() . the condition 2z{j + k) =7 



gives 



1-y 



2[1- {l + 4TrK^)y] ' 



(3.14) 



Since j + k = k/(1 — 2z), this corresponds to the positive 
(negative) branch of Eq. H3.11|l for z < (>)l/2. These 
three conditions for the gradient of t he sc alar field to be 
null are shown in Fig. O for n = l/^/T2n. Note that the 
positive and negative branches of Eq. H3.11|l should not 
really be represented in the same yz diagram, since the 
conditions 7 = and 2z(7 -I- k) = 7 clearly cannot be 
satisfied simultaneously, but we do so for convenience. 



If the gradient is spacelike, the energy density e given 
by Eq. H2.3|l is negative. This is not unphysical (eg. the 
weak energy condition is always satisfied), even though 
the identification with a stiff fluid is no longer possible. 
However, 7 must be real and, from Eq. H3.10|) . this applies 
within three regions in the yz plane: 



or 



or 



z < — and 



1 



y> 



1 



1 + Attk^ ' 



z < — and y < 0- 



These are the three unshaded regions in Fig. O Ahhough 
one might speculate that the complex 7 region could be 
associated with a complex scalar field, this cannot be 
confirmed within the context of the present single-field 
analysis. 



B. Similarity and trapping horizons 

From Eqs. H3.1|l and (|3.2|l . the induced metric on the 
surface ^ = const is 



ds^ 



-y-^z-^{l-2z)dr^+r^{d0^+sm^9d(t>'^). (3.15) 



Thus, providing y y^ 0, z = 1/2 corresponds to a simi- 
larity horizon and this is either a particle horizon or an 
event horizon. Although z = 1/2 plays a similar role in 
the field description to T^^ = 1 in the fluid description, 
we will show that the degeneracy in the fluid description 
is resolved in the field description. From Eq. 13.10|) . 7 
can be finite at z = 1/2 only if y = 1/(1 -I- Attk'^). 

If y = 0, one has a trapping horizon. This notion 
was defined in [33, llg and has turned out to be a very 
useful generalization of the concept of apparent hori- 
zons. Past trapping horizons include cosmological ap- 
parent horizons, while future ones include black hole ap- 
parent horizons. Equation (|3.1U|I implies that either z or 
7 must diverge when j/ = 0. From Eqs. 1)3. l|l and H3.2|l . 
dr/dv = g/2 = e^/(2z) is an outgoing null surface, so a 
future trapping horizon has infinite \z\ and finite 7, while 
a past trapping horizon has z — 1/2 and infinite 7. More 
precisely, y oc z~^ in the vicinity of future trapping hori- 
zons, while y ex (1 — 2z) and 7 oc y~^ in the vicinity of 
past trapping horizons, where we have used Eq. I|3.10|l . 
In the latter case, the apparent singularity arises from 
the fact that ^ or r has an extremum along ingoing null 
rays v = const on past trapping horizons. 

From these properties of similarity and trapping hori- 
zons, we can deduce the following. A similarity hori- 
zon with z = 1/2 cannot be a future trapping horizon. 
z — 1/2 is satisfied at both similarity and past trapping 
horizons. We can show that a black hole event horizon is 
untrapped with negative ingoing null expansion and pos- 
itive outgoing null expansion, where we need to assume 
the existence of future null infinity for the definition of a 
black hole event horizon. This is because, if both expan- 
sions were negative, no causal curve in the past neigh- 
bourhood of the black hole event horizon could reach 
future null infinity. On the other hand, if both expan- 
sions were positive, there would exist a causal curve in 
the future neighbourhood of the black hole event horizon 
which could reach future null infinity. It also follows that 
a black hole event horizon cannot coincide with a past 
trapping horizon. Therefore, a black hole event horizon 
is untrapped, i.e. its null expansions have different signs. 



C. Transformation between the two formulations 

The key step is to obtain the correspondence between 
two equivalent expressions for the three quantities 2m/r, 
Sner"^ and V. The first two can be written as 

2m/r = AEIVI"^ ^l~y, (3.16) 

Sner^ ^ ^1^1"^ = 47r7y[2z(7 + k) - 7]. (3.17) 

To obtain the expression for V in the field description, 
we introduce the following vector fields: 



\a, "a 



^a^,b, 



(3.18) 



where eat is the totally antisymmetric tensor, which sat- 
isfies e^ecb — —gab- Since these two vectors satisfy the 
normalisation relation 

Mqu" = — fifcfi , (3.19) 

the velocity function V can be written as 



V 



U'^ia 



(3.20) 



if (/i.a is timelike; here ^° = (9^)°. Since ^ is a function of 
X, this definition applies for any self-similar spacetime 
with a scalar field. Calculating the right-hand side of this 
equation in the Bondi coordinates and assuming k 7^ 0, 
we have 



Z I K 



(3.21) 



where e^^ is chosen to be positive. These relations give 
the explicit coordinate transformations between the fluid 
and field descriptions if ip^a is timelike. 

Equation (|3.21() means that the condition V^ — 1 in 
the fluid description separates into two independent con- 
ditions, z = 1/2 and 7 = 0, in the field description. 
From Eq. 13.15|) . the condition z — 1/2 implies that the 
^ = const surface is null unless y — 0- On the other hand, 
from Eq. (|3.12|) . the condition 7 = implies 0.a is null 
and this is plotted in Fig. fusing Eq. H3.13|) . There is no 
coordinate singularity even when (jj^a is null in the field 
description, and this is very important for the present 
analysis. As 7 approaches zero, Eqs. 1)2. 811 . 13.17|l and 
(|3.21l) imply that 1 — V^^ is proportional to 7 and that 
e2''(l _ v^) has a finite limit e^ ^(l - 2z)/{yz'^). This 
is zero only if z = 1/2, so Eq. H2.14|l implies that the 
similarity surface is null only in this case. 



D. Friedmann solution 

We now consider the flat Friedmann solution. This has 
v4 = 1/3 everywhere in the fluid description and a regular 
simflarity horizon. Since Attk'^ = 1/3, Eqs. ^J^~^J^ 
imply 



y2^^2^4X/3 ^ 



E ■ 



:C'e 



3„2X 



(3.22) 



The integration constant C conies from the remaining 
rescahng freedom of t. Using Eqs. l|TTU|) . (|TTO|) . l|TT7jl . 
H3.21(l and H3.22(l . it can be shown that the Friedmann 
solution is described by: 



^ n _L o /i r ~ expansion , 3.23 

2/(1 + 2^1-2/) 



2/(1 + 2VT^) 



(collapse). (3.24) 



The discussion about trapping horizons confirms that 
the first and second solutions represent the expanding 
and collapsing Friedmann solutions, respectively. The 
expanding solution is shown by the thick solid line in 
Fig. 121 and this crosses z = 1/2 at y = and y = 3/4. 
From Eq. (|3.21l) V ^ 1 a,t y = 3/4 and so this corre- 
sponds to the cosmological particle horizon. Eq. (|3.23|l 
implies y~^{l — 2z) -^ —1/12 as y ^ 0, which leads to 
V — 2 from Eqs. H3.11|l and H3.21|l . Therefore this surface 
is spacelike and corresponds to a past trapping horizon. 
This is the cosmological apparent horizon, which is out- 
side the particle horizon and coincides with the Hubble 
horizon Q. This illustrates that the condition z — 1/2, 
like the condition V = 1, does not always correspond to 
a null surface. The causal structure of the Friedmann 
solution in this case is shown in Appendix lEl 



E. Structure of similarity horizons 

As we have seen, a similarity horizon is defined as a 
null similarity surface. It is denoted by ^ = ^s and char- 
acterized by z = 1/2 and y ^ 0. Since this corresponds 
to a singular point of the autonomous system (|3.7|l - l|3.9|l . 
the usual uniqueness of solutions may break down there. 
The behaviour of solutions around ^ = £s can be anal- 
ysed by dynamical systems theory 16^, 32, 39J. First, 
we restrict our attention to k 7^ and assume that 7 
has a finite limit at ^ = ^s- Then (1 — 2z)'j tends to 
zero at the similarity horizon from Eq. (|3.9|) . Otherwise, 
7 oc In 1^ — ^s| in the vicinity of ^ = ^s and this necessarily 
implies a spacetime singularity unless y = 0. 

The full details of the analysis are described in Ap- 
pendix ^ Here we describe the important qualitative 
results. The regularity of the similarity horizon requires 
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Vs = 



1 + Attk'^ 



1 



(3.25) 



from Eqs. igSJ and l|TTn|l . We linearise the ODEs 
around the similarity horizon and represent the solutions 
as trajectories in (z,y,j,^) space. Generically regular 
solutions can only cross the similarity surface along two 
directions in this space, corresponding to two eigenvec- 
tors ei and 62. The associated values of 7s are: 
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(1 -a){l + a) 
' a^(l-2a) ' 



72 = ±00 



(3.26) 



where a = A-kk? . Along the second eigenvector, there is 
a solution which describes a density discontinuity surface 
at fixed ^ = ^s- We call this a shock- wave solution. 

As discussed in Appendix^ the signs of the eigenval- 
ues associated with the two eigenvectors determine the 
qualitative behaviour of solutions around the equilibrium 
point. For < AtthP' < 1/2, the similarity horizon is a 
non-degenerate node, with ei and 82 corresponding to 
the secondary and primary directions, respectively. This 
means that there is a one-parameter family of solutions 
which belong to the 62 direction and an isolated solu- 
tion which belongs to the ei direction. In particular, 
this includes the case with Attk^ ~ 1/3. Since the flat 
Friedmann solution is along ei, this is isolated, while all 
other solutions have a sound-wave with diverging 7 at 
the similarity horizon. One of them is a shock-wave. For 
Attk'^ = 1/2, the similarity horizon is a degenerate node. 
For 1/2 < 4ttk^ < 1, it again becomes a non-degenerate 
node but with ei and 82 corresponding to the primary 
and secondary directions, respectively. For Attk'^ > 1, 
the similarity horizon is a saddle point. This means that 
there are two isolated solutions associated with ei and 

62. 

Brady ^j| has analysed the structure of the system 
governed by Eqs. (|3.7|I - H3.1U|I for all values of the con- 
stant K. For < Attk^ < 1, which includes our choice 
of Attk^ = 1/3, the allowed region in the yz plane is 
similar to that shown in Fig. The solution can cross 
z = 1/2 with finite 7 only if ys and 7s are given by 
Eq. H3.25(l . Brady has also made a non- linear analysis 
of the similarity horizon and we expand his analysis in 
Appendix ^ correcting a typographical error. A com- 
parison of the two analyses shows that the qualitative 
behaviour obtained by linearizing the ODEs reflects the 
true non-linear behaviour. 



F. Non-existence of self-similar PBHs with a 
massless scalar field 

If the initial perturbation extends to infinity, so that 
the specific binding energy is negative there, the pertur- 
bation never goes to zero because the angular part of the 
metric acquires a numerical factor which results in a solid 
angle deficit |2a|. Although these solutions are some- 
times described as asymptotically Friedmann 3; they 
could in principle be observationally distinguished from 
the exact Friedmann solution by the angular diameter 
test. However, inflation itself would not produce models 
of this kind. We therefore first seek a black hole solution 
surrounded by an exact flat Friedmann solution beyond 
the sonic point, which is also the particle horizon, and 
consider matching to other interior self-similar solutions 
there. 

The required behaviour of the function z for a black 
hole is shown by the solid line in Fig. 0] As we move 
along the ingoing null surface v = const from the cos- 
mological apparent horizon, ^ decreases from ^cah- Thus 



z(^) increases from 1/2, reaches a maximum, turns to 
decrease and then crosses 1/2 at fph- Up to this point, 
the solution is described by the Friedmann solution, but 
the uniqueness of the solution breaks down at ^ph. Be- 
cause the Friedmann solution approaches the similarity 
horizon along a secondary direction, it has to connect to 
a member of the one-parameter family of solutions as- 
sociated with the primary direction. As we decrease ^ 
further, in order to have a black hole solution, z should 
first decrease, reach a minimum, begin to increase and 
then cross 1/2 at the black hole event horizon ^ = ^beh- 
It should then diverge at the black hole apparent hori- 
zon ^ = ^ah < Cboh (which is necessarily within the event 
horizon), although this feature is not required for the 
following discussion. 

It is easy to see that the required behaviour is impos- 
sible. Because z < 1/2 is required for ^bch < S, < Cph, 
we need y > 3/4 to avoid the shaded region in Fig. 13 
and Eq. H3.7|l then implies 2/3 < z/z < 2. However, this 
means that z cannot have an extremum, which gives a 
contradiction. The solution therefore either enters the 
negative mass region beyond y = 1 or approaches z — 
as S. —* — oo (i.e. as r — > 0), corresponding to the ex- 
act Friedmann solution. The only way to invalidate this 
proof would be to find a solution which goes through the 
lower shaded region in Fig. 13 One can never reach the 
black hole event horizon by going around this region. 

Figure |2| shows solutions obtained through numerical 
integration of Eqs. H3.7|) - H3.9|l inside the particle hori- 
zon. Some solutions directly cross y = 1 and go into the 
negative mass region with y > 1. However, there are 
also solutions which reach y = 3/4 at z < 1/2. These 
necessarily cross the curve 7 = 0, which means that they 
have the required behaviour for V'^ in Fig. [5] but not for 
z in Fig. 21 One could also consider solutions which pass 
smoothly into the region where the scalar field gradient 
is spacelike. However, in this case, the above proof still 
applies because z remains monotonic, so there is again 
no black hole solution. Such solutions have to touch the 
line y = 3/4, after which they become negative-branch 
solutions and eventually enter the negative mass region 
with y > 1. 

Although we have assumed that the spacetime exterior 
to the black hole solution is an exact Friedmann solution, 
the proof also applies if it is asymptotically Friedmann. 
For the equality Attk^ = 1/3 must still be satisfied in 
order to allow the exact Friedmann solution, the black 
hole solution must still have a particle horizon at z = 1/2 
and there must still be a region where < z < 1/2 inside 
the particle horizon. Therefore one can still conclude 
that z cannot cross z ~ 1/2 again from below. This 
proves that there is no black hole event horizon inside 
the particle horizon even for asymptotically Friedmann 
solutions. 

It should be noted that the situation is quite different 
for a perfect fluid with equation of state p = kp (0 < k < 
1). In fact, in this case, analytic and numerical studies 
show that there do exist asymptotically Friedmann self- 
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FIG. 4: The required form of 2;(^) for a black hole solution 
is plotted as a solid line if the cosmological apparent horizon 
is spacelike and as a dashed line if it is timelike. 



similar solutions with a particle horizon and black hole 
event horizon 0, [la, uM ■ The crucial point is that there 
is no extremum in z as a function of ^ in a subsonic 
region but there can be in a supersonic but subluminal 
(z < 1/2) region for < fc < 1. The reason why the 
non-existence proof works in the scalar case is that the 
sonic horizon in the scalar field system coincides with the 
similarity horizon. Since z cannot have an extremum in a 
subluminal region, the non-existence of self-similar black 
holes in the scalar field system follows immediately. 

If there is a shock-wave, it can be located only at the 
similarity horizon. As we have seen, such a solution has 
z = 1/2 and y = (1 -I- 47r/c^)~-^ but changing 7. If we 
consider a solution which is interior to this shock-wave 
and assume a regular junction, it must also have z — 
1/2 and y = (1 -I- 47rK^)~"'^ there. Because of regularity, 
7 must have the same value (47rK)~^ at the similarity 
horizon and in the interior solution. In fact, this shows 
that no regular shock-wave is allowed. We can still have a 
sound-wave solution with diverging 7 for < Attk^ < 1/2, 
which is included in the above proof. If we do not assume 
regularity of the similarity horizon, we may consider an 
extension beyond the similarity horizon to a solution with 
a different k. However, the different value of k means a 
different value of y for z = 1/2 because y = (1 -I- 47rK^)~^ 
there. One needs a singular null hypersurface because of 
the discontinuity in the Misner-Sharp mass. We do not 
consider this possibility further here. 



G. Bicknell and Henriksen's PBH construction 
with a null dust 

Bicknell and Henriksen [23 found an extension of the 
self-similar stiff fiuid solution beyond the timelike V^ — 1 
surface, containing null dust. The extension is given by 



an ingoing Vaidya solution 40] , in which the hne element 
is given by 



ds' 



1 



2ni(v) 



jdv^ + 2dvdr + r^{dd^ + sin^ 6'#2 ) . 



(3.27) 

The stress-energy tensor for null dust can be written in 
the form 



Tab — Crlalb, 



(3.28) 



where la is a null vector. For the Vaidya solution, la — 
—(5!? and a is written as 



1 dm 

47rr-^ dv 

The similarity assumption implies 

m{v) — /if, 



(3.29) 



(3.30) 



where [i, is an arbitrary constant. Then, in terms of y, z 
and ^, Eqs. (|3.2|) and (|3.1t)|) imply that this solution is 



y 



Eliminating ^, we have 



1 - 2/ie-«, 



1 - 2/ie-« ' 
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(3.31) 
(3.32) 

(3.33) 



The junction between the Vaidya metric and the stiff fluid 
solution can be implemented on the line 7 = 0, given by 
Eq. H3.13|l . For 47rK^ = 1/3, we then have 



/^ 



-(1-y*)', 



(3.34) 



where y = y* on the junction surface. Thus < /i < 3/64 
for 3/4 < y, < 1. We can see explicitly from Eqs. (|3.7() 
and H3.8|l that z and y are continuous on this surface, 
which guarantees the continuity of the first and second 
fundamental forms of the junction surface. For the ex- 
tended solution to have a black hole event horizon, we re- 
quire that z go below 1/2 for y > and Eq. (|3.33|) shows 
that this condition is necessarily satisfied. An example 
of an extended solution is shown by a thin dashed-dotted 
line in Fig. |31 

If the massless scalar field </> has a null gradient, the 
stress-energy tensor (|2.2|l is equivalent to that for the null 
dust (|3.28() . with la = a^^l'^XI a<i>- For this reason, one 
might expect that this construction could be realised con- 
sistently with a massless scalar field alone. Nevertheless, 
this is not true. If a null dust is equivalent to a massless 
scalar field, the vorticity-free condition must be satisfied: 

V[a(\/^^6])=0. (3.35) 

However, for the Vaidya solution, we can see 



V[o(\/ct?i] 



(3.36) 



This means that the vorticity-free condition is not satis- 
fied for the Vaidya solution unless it is vacuum. There- 
fore, the null dust which appears in the Vaidya metric 
cannot be realized by a massless scalar field. 

The stiff fluid or massless scalar field in this construc- 
tion turns into null dust on the timelikc surface where 
7 = 0, but it is then no longer describable as a massless 
scalar field. The Bicknell and Henriksen PBH solution 
exploits this feature and that is why this solution goes 
through the prohibited region where 7 is complex and 
why it circumvents the above proof that z cannot have a 
minimum. However, this may be regarded as physically 
contrived. 



IV. SCALAR FIELD WITH A POTENTIAL 

A. Autonomous system 

When the scalar field has a potential ^((/'), the ex- 
istence of self-similar solutions requires it to be of expo- 
nential form [4ll |: 

y(0) = V^e'^^'^. (4.1) 

Then the stress-energy tensor of the scalar field is 

'1 



r„, 



h- gah\ 7:4>,c4>''' + V^e 



SirXcti 



(4.2) 



With the self-similarity ansatz, the Einstein equations 
and the equations of motion for the scalar field reduce to 
the following ODEs: 

(Ing)- = inh , (4.3) 



9 



g[l~87Tx'Voe^^''j--g, (4.4) 
8n(h + k) [x(h + k) - gti] , (4.5) 



(4.6) 



{ghy+gh = VS^Xgx^Voe'^^'^^'' 
+2x (h + h + n) , 

where we use the same notations as for the massless case, 
except that one requires k — (y/ZTrX)^^ in the presence 
of a potential. This system has been investigated in the 
context of the kink instability analysis [42| . 

We define the functions y, z and 7 in the same way as 
before [cf. Eq. H3.2|) ] and introduce a new function 



P = SttVo exp ( 2^ H 

Eqs. (|3.7f) - (|3.10|) are then replaced by 

i = z[2-y-\l~P)], 
V = l-/3-(4^72 + l)y, 

(l-2z)7 = 2Kz-7[y"i(l-/3)-2z] 



/? 



P 



2^(^+3 



47rKy 



(4.7) 



(4.8) 
(4.9) 

{^.10) 



(4.11) 
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with the constraint equation 

y[{l + Attk^) - 47r(7 + Kf{l - 2z)] = 1-/5. (4.12) 

The sohitions are therefore constrained to he on a 3- 
dimensional hypersurface in a 4-diniensional phase space. 
Because the dimension of the phase space is higher than 
that for the massless case, the analysis becomes more 
comphcatcd. However, the number of degrees of freedom 
for the self-similar solutions is the same because the pa- 
rameter K is fixed by the potential, whereas it is freely 
chosen for the massless case. 

B. Friedmann solution 

The flat Friedmann solution satisfies these equations 
providing k = {'\/2t:X)~^. This is in contrast to the 
massless case, in which the Friedmann solution requires 
Attk^ = 1/3. In the presence of a potential, the scale 
factor is given by 



i{t) ex t^ 



ex t 



2/A2 



(4.13) 



Thus the expansion is decelerated for < Atth"^ < 1, 
accelerated for Attk^ > 1 and with a constant velocity 
for Attk^ = 1. We call the scalar field in the accelerated 
case a "quintessence" field. As discussed in Appendix IbI 
there are five types of causal structures, corresponding to 



< 47rK2 < 1/2, 471^2 
and 4:'kk'^ > 1. 



1/2, 1/2 < Attk^ < 1, Attk^ = 1 



C. Structure of similarity horizons 



As seen in Section IlII Bl similarity horizons are char- 
acterised by z = 1/2 and y 7^ 0, future trapping horizons 
by y = and \z\ = 00, and past trapping horizons by 
2/ = and z = 1/2. From Eqs. (|^?Tnjl and ^J^, the 
regularity of the similarity horizon yields 



Vs 



1-ft 
1 + Attk^ 



j3s + iTTn'^ 



(4.14) 



where we have assumed /3s 7^ 1- Eq. 14.1211 shows that 
/3s = 1 is only possible when 7 = 00, which means that 
the similarity horizon coincides with the past trapping 
horizon. Note that Eq. (|4.14() reduces to Eq. H3.25|l when 
/3s -0. 

As is before, we linearise the ODEs around the similar- 
ity horizon, which again corresponds to a singular point 
in the presence of a potential. The situation is qualita- 
tively similar to the massless case, except that the solu- 
tions are now represented as trajectories in (z, j/,7,^,/^) 
space. Generically regular solutions can only cross the 
similarity surface along two eigenvectors ei and 63. As 
before, these are associated with a finite and infinite 
value of 7s, respectively. We do not give the analogue 
of Eq. (|3.2()|) explicitly, since it is very complicated, but 



the infinite value is associated with a shock-wave solu- 
tion, for which y ~ Hs, z = 1/2, /3 = /3s, ^ = ^s but 
7 changes. This belongs to the one-parameter family 
of solutions which cross the similarity horizon along 62. 
As before, the similarity horizon behaves as a node for 
< AiTK?' < 1 or as a saddle for Attk? > 1. Appendix IFl 
provides the details of the linearised ODEs and eigenvec- 
tors for this case. 



D. Non-existence of self-similar PBH in 
decelerating Friedmann universe 



We now consider solutions which are exactly Fried- 
mann outside the particle horizon. The Friedmann back- 
ground is decelerating for < Attk'^ < 1, i.e. for a steep 
potential with A^ > 2. As shown in Appendix IbI the par- 
ticle horizon in this case is a similarity horizon and the 
Penrose diagram for a solution containing a black hole 
would be as indicated in Fig. ^ Because of the nodal 
nature of the similarity horizon, there is a one-parameter 
family of interior solutions which can be matched to this 
background. However, we now show that no interior self- 
similar solution can contain a black hole event horizon. 

For < Attk^ < 1/2, i.e. for A^ > 4, the cosmologi- 
cal apparent horizon is spacelike and outside the parti- 
cle horizon, as shown in Appendix ^ So if we have a 
black hole event horizon, as ^ decreases monotonically 
along a future-directed ingoing null ray from the particle 
horizon, z decreases from 1/2, reaches a minimum, in- 
creases and then crosses 1/2 again, y is positive in this 
region. The required behaviour of z is therefore as plot- 
ted with a sohd curve in Fig. ^ For 1/2 < inn'^ < 1, 
i.e. 2 < A^ < 4, the cosmological apparent horizon is 
timelike and inside the particle horizon, as shown in Ap- 
pendix^ Since a black hole event horizon is untrapped, 
as ^ decreases monotonically along a future-directed in- 
going null ray from the cosmological apparent horizon, 
z decreases from 1/2, reaches a minimum, increases and 
then crosses 1/2 again, y increases from zero as we move 
along a future-directed ingoing null ray from the cosmo- 
logical apparent horizon. The required behaviour of z is 
plotted with a dashed curve in Fig. 21 For Attk^ = 1/2, 
i.e. A^ — 2, the cosmological apparent horizon coincides 
with the particle horizon, as shown in Appendix IBI 

We now show that such a solution cannot exist. From 
Eqs. (|4.8|) and 14.12|l . we have 



= 1 - Attk^ + 47r(7 + k)2(1 - 2z). 



(4.15) 



Therefore, if Attk'^ < \, z cannot have an extremum for 
z < 1/2. Hence there is no self-similar black hole solution 
in this case. The proof also applies for the asymptotic 
Friedmann background. Thus the conclusion that there 
is no self-similar black hole solution in the case without 
a potential also applies for a decelerating universe. 
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black hole singularity • 




FIG. 5: The causal structure of self-similar black holes in 
an accelerating universe. The trajectories of similarity sur- 
faces are shown. There is a null big bang singularity, a null 
naked singularity and a spacelike black hole singularity. The 
cosmological event horizon (CEH) and the black hole event 
horizon (BEH) are both null similarity surfaces, i.e., similarity 
horizons. 



E. Non-existence of self-similar PBH in 
accelerating universe 

If Attk'^ = 1, or A^ = 2, the Friedmann universe ex- 
pands with a constant speed. In this case, the big bang 
singularity is an outgoing null surface and future null in- 
finity is ingoing null surface, as shown in Appendix^ We 
proved in Appendix IbI that there is no similarity horizon 
in this case, so we cannot match the exterior Friedmann 
background to an interior black hole solution. We there- 
fore exclude this case. 

For inn'^ > 1, i.e. for a flat potential with < A^ < 2, 
the scalar field with a potential acts as dark energy, which 
causes accelerated expansion. The Penrose diagram for 
a self-similar black hole spacetime is now different from 
the decelerating case, as indicated in Fig.[Sl This reflects 
the fact that the causal structure of the flat Friedmann 
universe is itself modified, as discussed in Appendices El 
and[ni The spacetime now has a null big bang singular- 
ity, as well as a null naked singularity, a spacelike black 
hole singularity and black hole event horizon inside the 
cosmological event horizon. The similarity surfaces are 
spacelike outside the cosmological event horizon and in- 
side the black hole event horizon, timelike in between 
them and null on the horizons themselves. So the sim- 
ilarity horizons can be identified with the cosmological 
event horizon or black hole event horizon. 

In this case, the sign of the right-hand side of Eq. (|4.15|) 
is apparently indefinite, so we cannot immediately apply 
the earlier proof. However, the nature of the similarity 
horizon in the Friedmann solution now plays a crucial 
role. Because the cosmological apparent horizon is time- 
like and inside the cosmological event horizon, as shown 
in Appendix ^ the external Friedmann solution could 
only be matched to an internal self-similar solution at the 
cosmological event horizon. As shown in Appendix ^ 



this horizon is a saddle, so there are just two isolated 
solutions which can be matched to the Friedmann back- 
ground solution: one is Friedmann itself and the other is 
a shock-wave solution. However, since we assume regu- 
larity of the cosmological event horizon, z, y and /? must 
be the same in the shock-wave solution and the exterior 
Friedmann background. The regularity therefore requires 
the same value of 7 and hence this similarity horizon 
is identical to that of the exterior Friedmann solution. 
Since the interior can only be the Friedmann solution it- 
self, we can conclude that there is no self-similar black 
hole solution interior to the Friedmann background. This 
proof does not exclude the possibility of a self-similar 
black hole in the asymptotically inflationary Friedmann 
back gro und and this possibility is now under investiga- 
tion 



V. CONCLUSION 

In summary, we have shown that there is no self-similar 
black hole solution surrounded by an exact flat Fried- 
mann solution for a massless scalar field or a scalar field 
with a potential. This extends the previous result ob- 
tained for a perfect fluid with equation of state p = ke 
(0 < fc < !)• Indeed we can summarize this important 
result as a formal theorem. 

Theorem. Let (M,g) be a spacetime which satisfies the 
following four conditions: (i) it is spherically symmet- 
ric; (a) it admits a homothetic Killing vector; (Hi) it 
contains a stiff fluid or a scalar field with non-negative 
potential; (iv) it contains no singular hypersurface. Let 
{A4,g) also satisfy at least one of the following two con- 
ditions: (v) it coincides with the flat Friedmann solution 
outside some finite radius; (vi) it is asymptotic to the 
decelerating Friedmann solution. Then (M,g) has no 
black-hole event horizon. 

This result disproves recent claims that PBHs could 
grow self-similarly through accreting a quintessence field. 
It also suggests that accretion onto black holes is sup- 
pressed by the cosmological expansion in a scalar field 
system. This is consistent with numerical simulations, 
which show that PBHs with different initial sizes even- 
tuall y b ecome much smaller than the particle hori- 
zon [23, y^. Note that, for the massless case, it has 
already been shown analytically that the accretion rate 
is exactly zero for a black hole whose size is exactly that 
of the cosmological apparent horizon _29, _30j . 

Although Bicknell and Henriksen's construction shows 
that there are self-similar black hole solutions satisfy- 
ing the necessary energy conditions, such solutions are 
ad hoc and probably not physically realistic. Certainly 
the PBH mass could grow self-similarly (M ex t) only 
with very special matter fields. On the other hand, if we 
drop self-similarity condition, there might still be black 
hole solutions which exhibit an interesting amount of 
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accretion. For example, kinematically self-similar solu- 
tions |43 . |45| contain characteristic scales and this sug- 
gests consideration of a wider class of matter fields. In 
this case, one might have solutions in which the PBH 
grows as M oc i^'" {a ^ 1) for some similarity index a. 
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APPENDIX A: COMPARISON WITH OTHER 
WORK 

It is interesting to explain why the results of this paper 
disagree with those presented in various earlier papers. 
First, we clarify the flaw in the claim by Lin et al. |2g 
that a self-similar black hole solution can be attached to 
the Friedmann particle horizon if the equation of state 
is stiff. This flaw was first pointed out by BickncU and 
Henriksen [23 but our analysis has illuminated its na- 
ture. Lin et al. use different variables from those in 
Section ^] but it is easy to express the results in our no- 
tation. They prove analytically that, as the similarity 
variable X decreases, the function V must reach a mini- 
mum below 1 and then rise again to 1 (cf . Fig. ^ . They 
correctly point out that the pressure and velocity gradi- 
ent diverge at this point but they fail to appreciate that 
the density function E goes to zero and that the metric 
function i' goes to infinity, so that the surface V^ — 1 
is timclikc rather than null. They therefore misinterpret 
this point as a black hole event horizon. They argue that 
the divergence of the pressure gradient is non-physical, 
since it can be removed by introducing an Eddington- 
Finkelstein-type coordinate, but do not notice that the 
pressure gradient also diverges at the cosmological par- 
ticle horizon. However, it is worth stressing that Lin et 
al. do correctly deduce the form of the solution up to the 
inner point where V^ = 1. 

Second, we point out the source of our disagreement 
with the claim of Bean and Magueijo H| that black 
holes can accrete quintessence fast enough to grow self- 
similarly. They consider a scalar field in a potential of 
the form given by Eq. I|4.1(l . By generalizing the analysis 
of Jacobson |4g, which attaches a Schwarzschild solu- 
tion to a cosmological background in which the field has 
the asymptotic value 0c, they claim that the energy flux 



through the event horizon is T^y = 
hole accretion rate 



dM 

~dr 



16ttM^ 



leading to a black 



(Al) 



They correctly infer that only the kinetic energy of the 
scalar field is accreted, the scalar potential merely influ- 
encing the evolution of the asymptotic background field: 



AV87r 



log 



(A2) 



where to is a constant of integration. The accretion rate 
therefore becomes 



dM KAP 



dt t^ 

with K = 8/A^, which can be integrated to give 

Mo 



M = 



^ _ KM^ /to 
to V t 



(A3) 



(A4) 



For Mq ^ K^^to, corresponding to an initial black hole 
mass much less than that of the particle horizon, this 
implies negligible accretion. However, for Mq = K~^to, 
Eq. I|A4|I predicts a self-similar solution, M ex f , in which 
the black hole always has the same size relative to the 
particle horizon. This is the result which Bean and 
Maguiejo exploit. For AIq > K^^to, the mass diverges at 
a time t = to/[l - to/iKMo)]. 

However, this analysis is exactly equivalent to that 
of Zeldovich and Novikov (3, which neglects the back- 
ground cosmological expansion and which Carr and 
Hawking disproved. Bean and Magueijo caution that the 
Carr-Hawking proof only applies for a perfect fluid and 
may not be relevant if the scalar field has a potential. 
However, the result of this paper confirms that there is 
no black hole self-similar solution even in this case. (Bet- 
twieser and Glatzel J4^ also question the Carr-Hawking 
analysis but they only use a general causality argument 
and do not give a detailed calculation.) It is interesting 
to note that a more refined Newtonian analysis (but still 
neglecting the background expansion) implies that the 
critical black hole size for self-similar growth exceeds the 
scale associated with a separate closed universe condition 
for an equation of state parameter exceeding 0.6 ;9]. This 
also indicates that the Newtonian prediction cannot be 
correct. 

Finally, we comment on the paper by Custodio and 
Horvath [3 . They also consider black hole accretion of a 
quintessence field with a scalar potential but, instead of 
Eq. HAljI. obtain 



dM 
~dt' 



27ttM^ 



(A5) 



where the numerical factor differs from that assumed in 
Eq. (jAip because of relativistic beaming. They disagree 
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with Bean and Magueijo's choice of the function 0c on the 
grounds that it neglects the local decrease in the back- 
ground scalar field resulting from the accretion. For rea- 
sons which are not altogether clear, they therefore prefer 
to focus on a model in which the quintessence flux into 
the black hole is constant. This requires that the poten- 
tial have a specific shape, which is not in fact of the form 
(14. 1|) required by self-similarity. This leads to the mass 
evolution 



M ^ 



Mo 



1 



KMo 
to 



(A6) 



(t-to) 



where X is a constant related to the (fixed) fiux. Again 
this is similar to the Zeldovich-Novikov formula. The 
mass diverges at a time t = to[l + l/{KMo)], just as 
Eq. (|A4|) does for Mq > K^^to, and they attribute this 
unphysical feature to the fact that the assumption of con- 
stant mass flux fails. They therefore consider alternative 
models in which the flux decreases as a power of time. 

Although we agree with their conclusion that the Bean- 
Magueijo argument is wrong, we would claim that Custo- 
dio and Horvath have not identified the more significant 
reason for its failure, namely that it neglects the back- 
ground cosmological expansion. In general, they find that 
the increase in the black hole mass is small unless its ini- 
tial value is finely tuned. However, according to the result 
in this paper, even the latter claim is misleading. 



APPENDIX B: GLOBAL STRUCTURE OF 
FRIEDMANN SOLUTION 

Since we are considering the possibility of black holes 
in a Friedmann background, we briefiy review the global 
structure of the flat Friedmann solution. This has the 
metric 

ds^ = -dt^ + a'^{t)[dx^ + x^{de^ + sin^ dcf)'^)] (Bl) 

with the scale factor a{t) = Ci^, where C and q are 
positive constants. The domain of the variables is < i < 
cxD and < X < oo- The big bang singularity is i = and 
the regular centre is x = 0. We have q — 2/[3(l-f'A:)] for a 
perfect fluid with p — kt {k > ~1) , q ^ 1/3 for a massless 
scalar field (fc = 1) and q = 2/A^ for a scalar field with 

the exponential potential V = V^e^^^^^'^ . The cosmic 
expansion is decelerated for < g < 1, accelerated for 
q > 1 and has constant speed for q = 1. 

In terms of the conformal time ry, defined by adrj = dt, 
we have a ex |7y|9/(i~9) for 9 7^ 1 and a (x e'^'' ior q = 1. 
The range of 77 is < 77 < 00 for < q < I, —00 < rj < 00 
for q = I and —00 < 77 < for q > 1. Defining new 
variables t and f by 



V 



tan 



t + f 



tan 



t + f 
tan I I — tan 



t-f 

2 
i-f 



(B2) 
(B3) 



we can draw the Penrose diagram for each case. Through 
this transformation, the big bang singularity t — and 
the regular centre x = are transformed respectively to 
i = and f = OforO<g<l and to i — f = — tt and 
f = for g > 1. It can be also shown lia] that future null 
infinity is given hy t + f — n for < q < 1 and t — for 
q > 1. The big bang singularity is spacelike for < g < 1 
but null for q > 1. 

When wc introduce the double null coordinates 



the line element takes the form 



V + X, 



ds^ = a^[-dudv + x^{dO^ + sin^ dd(j>^)], 



where 

^2 ^ ^2/(l-g) 

for q ^ I and 



(1-9) 



u + V 



2q/(l-q) 



^2^^2gC(n+.) 



(B4) 



(B5) 



(B6) 



(B7) 



for q — 1. The area radius r is given by r = a{v — u)/2. 
The Misner-Sharp mass m is written in this metric as 



2'^ 



^"^ ^U^ ^V 



Then 



2m 
r 



i-q 



for q y^ I and 



2to 
r 



C= 



{v-uY' 



(B8) 



(B9) 



(BIO) 



for q = I. The past trapping horizon, which coin- 
cides with the cosmological apparent horizon, is given 
by 2m /r = 1, or 



{2q^l)v 



for q ^ I and 



"-"^c' 



(Bll) 



(B12) 



for q — 1. Therefore, the cosmological apparent horizon 
is spacelike for < g < 1/2, null for q — 1/2 and timelike 
for q > 1/2. Figure El shows five different Penrose dia- 
grams, depending on the value of q. There is a particle 
horizon, given by u = 0, only for < q < 1. There is 
a cosmological event horizon, given by u = 0, only for 
q > 1. The flat Friedmann solution for q = 1 has neither 
a particle horizon nor a cosmological event horizon. 
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FIG. 6: The causal structure of the Friedmann universe with the scale factor a oc t'' for different values of q. The particle 
horizon (PH), cosmological apparent horizon (CAH) and cosmological event horizon (CEH) are indicated. 
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APPENDIX C: FRIEDMANN SOLUTION AS A 
SELF-SIMILAR SPACETIME 



APPENDIX D: ANALYSIS OF SIMILARITY 
HORIZON FOR THE MASSLESS CASE 



Although these flat Friedmann solutions are self- 
similar in the sense of homothety, it is not obvious in 
the present form. Using double-null coordinates, the 
standard form of spherically symmetric self-similar space- 
times is given by 

ds^ = -d.'^{a)dudd + uib^{a){dd^ +sm^ dd(l>^), (CI) 

where a = u/v is the similarity variable. Since the metric 
with the double null coordinates u and v does not satisfy 
this form, we transform to new double null coordinates 
u and V. 

For q 7^ 1, u and v are given by 

u = ±\u\^-^, v = ±\i\^^^, (C2) 

where the signs are chosen so that the sign of u (v) coin- 
cides with the sign of u (v) and we can show — 1 < cr < 1. 
Then the metric functions become 



d^ = const X 



const X 



-{l-q)/2 



±\<J 



{i-q}/2 



2g/(l-g) 



,(C3) 



|^|-(l-9)/2±|^|(l-«)/2 



2g/(l-g) 



x||ar(i-^)Tki^'"«f, 



(C4) 



where the upper and lower signs correspond to the posi- 
tive and negative values of a, respectively. For q = 1, 



and the metric functions become 



d — const, 



const X (Incr)^, 



(C5) 

(C6) 
(C7) 



where < cr < 1. 

The similarity surface is generated by a homothetic 
Killing vector. When this is null, it is called a similar- 
ity horizon and this does not depend on the choice of 
similarity variable. The induced metric on the similarity 
surface a ~ const is 



ds = —a (T du 



uib^{d0^ 



= -d^adv^ +uib^{d9^ 



sm 



sm 

2 



'^ (C8) 



'■)■ 



We provide two expressions because one of them may 
become invalid when one of the coordinates becomes de- 
generate. For g 7^ 1, the only similarity horizon is cr = 0. 
This is outgoing and null for < g < 1, corresponding to 
the particle horizon m = 0, and ingoing and null for q > 1, 
corresponding to the cosmological event horizon v — 0. 
For (7 = 1, there is no similarity horizon and the similar- 
ity surfaces are v — u = const, which are timelike. The 
similarity horizon is important because it is the matching 
surface between the external Friedmann background and 
the internal self-similar solution. 



We introduce a new independent variable u defined by 

(Dl) 



du 



where the similarity horizon ^ = ^s corresponds to u = 
±oo and is an equilibrium point. From Eqs. 13.7|l - (|3.9|) . 
we can write the ODEs near ^s in the form: 

J = z(l-2z)[2-y-i], (D2) 

^ = (l-2z)[l-(l + 4^7')y], (D3) 

du 

^ = 2Kz~j{y-^-2z), (D4) 

together with the constraint given by Eq. (|3.1U|) . We now 
write the asymptotic solution near ^g as 

z = -(H-xi), y = ys{l+X2), 7 = 73(1+2:3), ^ = Cs+a;4, 

(D5) 
where j/g and 7s are given by Eq. 13. 2511 . and regard 
(xi, X2, a;3, 0:4) as defining a 4- vector x. The linearised 
ODEs can be written in the matrix form 



du 



= Ax 



where 



/ a-1 

-1 
a — a 







1 -|- a 1 + a —a 

V -1 



(D6) 



(D7) 



and a = Attk^. This matrix generically has only two 
non-zero eigenvalues: 



Ai = —1 + a, A2 = —a. 
These are associated with two eigenvectors 



(D8) 



a(l -2a)(l-a) 

-(l-2a)(l-a)(l + a) 

{l-a){l + a) 

a{l-2a) 



62 



(D9) 

and the two associated values of 73 — {'^X'i/x4)s are given 
by Eq. (|3.2()|) . In particular, there is a solution belonging 
to the eigenvector 62 which has the form 



1 

^ = 2' y^y^" 



7 = 7s(l±e-n, e = 6, (DIO) 



and which we call a shock-wave solution. Note that the 
equilibrium point is a node for Ai < and A2 < 0, with 
the primary direction being associated with the higher 
value, and a saddle for Ai > and A2 < 0. 
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APPENDIX E: SIMILARITY HORIZON FOR 
THE MASSLESS CASE REVISITED 

Here we revisit Brady's non-linear anal ysis of the sim- 
ilarity horizon for a massless scalar field |31j . We intro- 
duce new variables ( and r? such that 



z^ ^ + C, y^Vs + v, 



1 



Vs 



1 + Attk'^ ' 



(El) 



where C — V = ^ corresponds to the similarity horizon 
and only ?]/( < is allowed. When ( and rj are suffi- 
ciently smah, Eqs. (|3.7|) - (|3.1U|) reduce to 



7 
drj 

dC 



-K± 



1 + Attk-^ 



8tt 



± 



2V8^ 



1 - Attk"^ 



1 

c 

1 -I- Atik^ rj 
1 - 47rK2 1 



(E2) 



(E3) 



Here and throughout Appendix El the upper and lower 
signs correspond to the positive and negative branches in 
Eq. I|E2|) . These equations can immediately be integrated 
to give 



± 



1 



27rK 



C\C\ 



i-KK^/il-A-ITK.^) 



(E4) 



7 = ± 

AlTK, 



1 , C(l -H 47rK2) 



ICI 



47rK^/(l-47rK^) 



where C is an integration constant. If ^ < and rj > 
0, Eq. IJE4|I with the upper sign coincides with Brady's 
solution (Eq. (C6) in 31j), except for the exponent of 
the second term on the right-hand side, which may be 
a typographical error. The negative- branch solutions do 
not cross the similarity horizon C = 77 = 0. 

The solution with the upper sign has the following be- 
haviour. If C = 0, it always crosses the similarity hori- 
zon. For < 47rK^ < 1, there is a one-parameter family 
of solutions which cross the similarity horizon, because 
the second term on the right-hand side of Eqs. I|E4|I and 
(|E5|) goes to zero as C ^ 0. For Ann"^ > 1, however, 
the solution with C = is isolated because 77 will not go 
to zero as C ^ if C 7^ 0. If we consider the deriva- 
tives along solutions with C 7^ 0, because of the second 
term, d-f/dC, is diverging for < Attk^ < 1/2 and fi- 
nite for 1/2 < Attk^ < 1. This reflects the behaviour of 
sound-wave solutions with C 7^ 0, which is discussed in 
Section Fill El using the linearised ODEs. Note that the 
shock- wave solution, which is given by Eq. (jDlOp . now 
has 77 = C = but changing 7. 



APPENDIX F: ANALYSIS OF SIMILARITY 
HORIZON FOR THE POTENTIAL CASE 

Using the variable u defined by Eq. Hm|l . the system 
of ODEs near the similarity horizon has the form: 



dz 
du 



z(l-2z)[2-y-i(l -/?)], 



(Fl) 



^ = (l-2z)[l-/3- (1 + 4^)7/], 

-— = 2kz 
du 



-i[y-Hl-p)-2z] + 



/? 



'iiTKy 



^ = 2(l-2z)/3(l-K7/^), 
du 

together with the constraint 

1 - /? = y[(l -I- Attk^) - 47r(7 + k)2(1 - 2z)]. 



(F2) 
(F3) 
(F4) 

(F5) 



The nature of the equilibrium point can again be classi- 
fied by the behaviour of solutions there. We put 

where xi to X5 are regarded as components of a 5-vector 
X. From the regularity condition at the similarity hori- 
zon associated with Eq. (|F3p and the constraint equation 
(|F5I) . we obtain Eq. 14.14|l . so that the similarity horizon 
is parametrized by the values of k and Ps ■ 

We now write the linearised ODEs in the form: 



du 



Ax, 



(F6) 



(Eoj where the matrix A is given by 



A = 



/ a-1 

A21 

-1 

V A51 



\ 
000' 

~a ^35 

000/ 



(F7) 



with 



A21 



^31 — ^32 — 
^35 



(1 + a)[l3s{l + a^) + a{l ~ a)] 

a^il-Ps^ 
x[{l-a)(3s+a], 



13s + a 
/3s{l + a)[(l-a)Ps + '2a] 



A51 



(l-/3s)(/3s + a) 

2(l + a)[(l-Q)/3, + a] 

a^l -f3s) 



(F8) 
(F9) 

(FIO) 

(Fll) 



where a = Attk^. This matrix has three zero eigenvalues 
and two generically non-zero ones, Ai = —1 -I- a and 
A2 = —a. Two of the eigenvectors associated with zero 
eigenvalues do not satisfy the constraint equation and 
so are ruled out. The eigenvectors associated with the 
eigenvalues Ai, A2 are respectively 



ei 



V (ei)^ / 



62 



1 





(F12) 



where 






{l + a)[{l + a^)f3s + a{l-a)] 





a3 


(l-/?s)2 




x[(l- 


a)/3s + 


a], 




(1 + 


«)[(!- 


a)Ps + a]Q 



a3(l-/3,)2(l-2a)(/3, + a)' 



(F13) 
(F14) 
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Q = -(l-a)(a2+a_i)/j2 

-2a{a^ + 2a-l)Ps + a'^il-a), (F15) 



5 _ 2(l + a)[(l-a)^s + a] 



(ei)^ = 



a^l-M 



(F16) 



These eigenvectors he on the constraint surface. 
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